Energy landscapes, ideal glasses, and their equation of state 
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Using the inherent structure formalism originally proposed by Stillinger and Weber [Phys. Rev. 
A 25, 978 (1982)], we generalize the thermodynamics of an energy landscape that has an ideal 
glass transition and derive the consequences for its equation of state. In doing so, we identify a 
separation of configurational and vibrational contributions to the pressure that corresponds with 
simulation studies performed in the inherent structure formalism. We develop an elementary model 
of landscapes appropriate to simple liquids which is based on the scaling properties of the soft-sphere 
potential complemented with a mean-field attraction. The resulting equation of state provides an 
accurate representation of simulation data for the Lennard- Jones fluid, suggesting the usefulness 
of a landscape-based formulation of supercooled liquid thermodynamics. Finally, we consider the 
implications of both the general theory and the model with respect to the so-called Sastry density 
and the ideal glass transition. Our analysis shows that a quantitative connection can be made 
between properties of the landscape and a simulation-determined Sastry density, and it emphasizes 
the distinction between an ideal glass transition and a Kauzmann equal-entropy condition. 



I. INTRODUCTION 

Deeply supercooled liquids and glasses occupy a promi- 
nent place in modern science and engineering. In addi- 
tion to their established and important presence among 
polymeric materials, glasses are becoming key elements 
in new technologies, such as pharmaceutical preservation 
and corrosion-resistant alloys jlj. The typical route to 
the vitreous state is to supercool a liquid fast enough so 
that crystallization is thwarted; hence the properties of 
glasses are intimately linked to those of the metastable 
liquids from which they are made At sufficiently low 
temperatures during the cooling process, the structural 
relaxation times of the metastable liquid become so slow 
that its mechanical properties begin to resemble those of 
a solid and the material is no longer in equilibrium rel- 
ative to a laboratory time scale. Though at this point 
its mechanical behavior is solid-like, the distinguishing 
structural feature of a glass is that it possesses no long- 
range microscopic order. 

The major industrial prominence of glasses tends to be- 
lie our rather limited theoretical understanding of these 
materials. Although significant progress has occurred 
in recent years [||, ^, ^, ^, many questions still re- 
main concerning the appropriate thermodynamic treat- 
ment of the glassy state [0, ^, ^, the relationship be- 
tween kinetics and thermodynamics in these systems 
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[TT[ |T|, [T|, |T|, lie), and the connection between 
molecular architecture and macroscopic behavior ^ pT[ . 
The present work addresses the first of these topics and 
aims to clarify some of the most conspicuous questions 
about the thermodynamics of supercooled liquids and 
glasses. The unresolved issues in glass thermodynamics 
are quite varied, and an important goal of this work is to 
investigate them with a common theoretical framework 
capable of complete thermodynamic description of the 
liquid state. We propose a simple landscape-based equa- 
tion of state for a supercooled liquid and demonstrate 
the usefulness of this approach by comparing theoretical 
predictions with simulation data. 

One unresolved issue regards the so-called Kauzmann 
temperature, Tk IQ |^, |18|, |1^]. As a liquid is supercooled, 
its larger heat capacity relative to the crystal causes the 
melting entropy to be gradually consumed until it ap- 
pears that at this distinguished temperature, the liquid 
and solid entropies become equal. Upon further cool- 
ing below the Kauzmann temperature, the liquid would 
eventually attain a negative entropy and would hence 
appear to violate the third law of thermodynamics. Ex- 
perimentally, however, the Kauzmann temperature must 
be extrapolated because upon cooling the glass transi- 
tion intervenes at higher temperature, thus preventing 
any such violation. Nevertheless for many liquids, the 
extrapolation needed to attain the equal-entropy condi- 
tion is quite modest, and attempts to understand the 
nature and implications of this impending entropy crisis 
underlie the thermodynamic interpretation of the glass 
transition. 

The idea that a kinetically-controlled glass transition 
prevents a thermodynamic catastrophe seems rather un- 
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settling and has led to the notion of a thermodynamic 
"ideal" glass transition at Tk (19|. In this sense, the 
experimentally observed glass transition is viewed as a 
kinetically-blurred manifestation of an underlying phase 
transition. The ideal glass transition has come to be as- 
sociated with the sudden entrapment of a system in a 
lowest-energy, unique amorphous configuration accessed 
in the limit of infinitely slow cooling. It is now recognized 
that in this definition, an ideal glass transition can ac- 
tually be a rather distinct occurrence from a Kauzmann 
equal-entropy point; the former is defined by the amor- 
phous state alone, whereas the latter makes reference to 
the crystal 

The energy landscape formalism introduced by Still- 
inger and Weber has been an important theoretical tool 
for formulating a thermodynamics of glasses . From 
this perspective, a liquid is described by the structurally 
distinct configurations through which it evolves, each 
termed an inherent structure, plus the kinetic "vibra- 
tional" distortions around these configurations. By defi- 
nition, the inherent structures are the complete collection 
of mechanically stable particle packings (local potential 
energy minima), and any one can be found from a given 
configuration by energy minimization. This permits a 
rigorous theoretical separation of liquid-state properties 
into inherent structure and vibrational components, in- 
cluding the energy, entropy, and pressure. The energy 
landscape paradigm has greatly facilitated the under- 
standing of low-temperature liquids and their glasses 

0, n, 0, 0, n, n, n, II, m. 

In the context of the energy landscape, ideal glasses 
acquire the rigorous definition of zero configurational en- 
tropy |l9) . The configurational entropy is that part of the 
entropy associated with the degeneracy of inherent struc- 
tures of a given potential energy; presumably when the 
temperature is low enough, the system samples only its 
minimum-energy amorphous conformation and its con- 
figurational entropy vanishes. In contrast, Kauzmann 
points are defined by the equality of total crystal and 
liquid entropies. Since the configurational entropy of a 
crystal is zero, the difference between an ideal glass tran- 
sition and a Kauzmann point is due solely to differences 
in vibrational entropies between the supercooled liquid 
and the stable crystal. In fact, it has been observed 
that several real substances do in fact exhibit Kauzmann 
equal-entropy points without violation of the third law, 
and in these systems the contribution of vibrational en- 
tropies is essential to the existence of a Kauzmann con- 
dition On the other hand, theoretical arguments 
have challenged the possibility of an ideal glass transi- 
tion at finite temperature by examining the effects of el- 
ementary excitations on the configurational entropy [ p^ . 
Still, whether an ideal glass exists in supercooled liquids 
as a thermodynamic phenomenon underlying the labora- 
tory (kinetic) glass transition remains an important open 
question. 

Our objective is to provide a simple landscape-based 
thermodynamic framework for liquids, including their 



equation of state, and apply it to the investigation of 
some of these questions. We make no attempt to de- 
scribe nonequilibrium states in our theory, such as the 
kinetically arrested glasses observed in experiment, but 
instead invoke the possibility of an ideal glass transition. 
In Section || we briefly review some of the empirical ob- 
servations made from computer simulations about the 
relationship between the equation of state and the en- 
ergy landscape. We proceed in Section III to outline our 
theory and examine some of its implications for these ob- 
servations as well as for the notion of an ideal glass tran- 
sition. We then describe in Section |^ what is perhaps 
the simplest quantitative model for an energy landscape 
of a simple liquid. Finally, in Section ^ we remark on the 
results that this equilibrium theory has for supercooled 
liquids and note areas of future research in light of them. 



II. THE EQUATION OF STATE IN THE 
LANDSCAPE PARADIGM 



Within the energy landscape formalism, it is possible 
to separate the equilibrium pressure of a liquid into con- 
tributions from its inherent structures and those due to 
vibrational displacements about the energy minima. The 
component of the pressure due to inherent structures. 
Pis, is particularly convenient to study because it can 
be measured directly by computer simulation. To cal- 
culate the inherent structure pressure at a given density 
and temperature, a large number of configurations are 
taken periodically from a molecular dynamics trajectory 
and their energy is minimized; this procedure locates the 
corresponding inherent structures, and the pressure of 
the minimized configurations is then calculated from the 
standard virial expression. Based on a number of in- 
vestigations which have used such a protocol. Pis has 
emerged as an important feature of supercooled liquids 
1^, 1^, m m, H, ||. In particular, it has suggested 
an unexpected connection between the liquid spinodal, 
microscopic heterogeneity, and the mechanical strength 
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The inherent structure pressure has been studied in 
computer simulations of several substances, and in all 
cases the density dependence of Pis has given rise to 
a picture qualitatively similar to that shown in Fig. |l| 
|l3|, |0|, P, H ||. It has been found in simple liq- 
uids that, above the melting line, the temperature of the 
equilibrium system from which the inherent structures 
are generated has very little effect on the type of curve 
shown in Fig. |l|, m eaning that Pis is approximately in- 
dependent of r [|13[. However, for all the systems studied 
so far below the melting temperature (and for more com- 
plex systems such as water even above the melting line! , 
there is a detectable influence of the temperature |2^, ^ . 

Fig. |l| reveals the presence of qualitatively distinct 
density regimes in Pis . As the density decreases from 
region III to II, the system's inherent structures evolve 
from a state of positive pressure to one of tension. The 



3 




FIG. 1: Schematic of the density dependence of the inherent 
structure pressure, Pis, for simple fluids. 



density at which the minimum in pressure occurs, the 
so-called Sastry density, ps , represents the point of max- 
imum attainable tension. At densities below this point 
(region I), the system fractures and its inherent struc- 
tures are no longer homogeneous but instead contain sig- 
nificant void regions . The Sastry density is therefore 
the density of limiting mechanical stability. In the case of 
glassy substances, whose mechanical properties are dom- 
inated by the slow relaxation from one inherent structure 
to another, the Sastry density offers an important mea- 
sure of their ultimate strength |l|, H, H, || ID . 

The study of the thermal and mechanical stability 
boundaries of supercooled liquids has received consider- 
able attention in recent years j2^, |4| ^ ^ . In a number 
of investigations, it has been found that the Sastry den- 
sity, determined from finite-temperature simulations, is 
nearly indistinguishable from the zero-temperature limit 
of the liquid branch of the spinodal curve, determined 
from accurate equations of state ||2^, U|. These results 
imply a very weak dependence of the Sastry density on 
the temperature of the equilibrated liquid from which it 
is calculated. Based on such results it has been proposed 
to interpret ps as the zero-temperature terminus of the 
liquid spinodal |Q. Under this hypothesis, the limit of 
mechanical stability of inherent structures coincides with 
the limit of thermodynamic stability of the liquid in the 
absence of thermal motion. This proposal suggests a par- 
ticularly striking connection between features of a liquid 
at high temperatures and absolute zero; however, the 
underlying ideas are controversial and they have yet 
to benefit from a rigorous formulation |2^, |34|. A par- 
ticularly prominent criticism of the spinodal connection 
stems from simulation results for a model glass-former in 



which the limit of mechanical stability for the liquid was 
explicitly calculated; in this counter-example the zero- 
temperature limit of liquid mechanical stability did not 
appear to converge to the Sastry density . 

This picture has become more intriguing with the re- 
cent observation in several simple models that the liq- 
uid spinodal and Kauzmann curve converge in the zero- 
temperature limit |2^, If the connection between 
the spinodal and Sastry density is made, this means that 
Tk{ps) = 0. The qualitative explanation for this oc- 
currence offered in |^ is that the inherent structures 
at Ps reach a point of zero configurational entropy, and 
hence represent an ideal glass, due to the constraint that 
they must be devoid of weak spots that would initiate 
fracture. That is to say, the requirements of homogene- 
ity and maximum tension limit the system to a single 
inherent structure at ps- This explanation invokes the 
notion of a system with zero configurational entropy; for 
this to apply to a Kauzmann point as well, the vibra- 
tional entropies of the crystal and metastable liquid must 
be equal. Though this is the likely behavior at abso- 
lute zero, more insightful theoretical or quantitative ar- 
guments substantiating the Kauzmann-spinodal connec- 
tion have not been offered. Furthermore, detailed calcu- 
lations for one particular model glassy system seem to 
exclude the zero-temperature spinodal-Kauzmann con- 
vergence 1 35 , questioning the universality of this obser- 
vation. 



III. THEORETICAL FRAMEWORK 

A. Free energy 

In the energy landscape formalism, the starting point 
in the analysis of a liquid is its potential energy surface 
(PES), the many-particle potential energy function of the 
configurational coordinates |Q . For a single-component 
system of N structureless particles, this surface exists 
in a space of 3N + 1 dimensions. In this enormous di- 
mensionality, the PES contains an overwhelmingly large 
number of local minima, the inherent structures, each 
corresponding to a mechanically stable packing of the 
particles |3^. By following a steepest descent trajec- 
tory, every configuration can be uniquely mapped onto 
one of these minima. Consequently, the entire config- 
uration space can be tiled into basins around inherent 
structures. The collection of configurations that corre- 
sponds to a given inherent structure is called its basin of 
attraction. This approach allows a rigorous separation of 
the partition function into contributions from inter-basin 
movement and vibration within a basin lE3l: 



-f3A 



/•0max(p) 

C / e 

"'0„,i„(rt 



7V[<T(0,p)-/30-/3avib(/3,0,p)], 



(3.1) 

where Z is the canonical partition function, (3 (= l/kT) 
is the Boltzmann factor, A is the Helmholtz free en- 
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ergy, C is a constant with units of reciprocal energy, 
4> is the potential energy per particle, p is the number 
density, cr((/), p) is the basin enumeration function, and 
o-\\h{Pi4>i p) is the basin vibrational free energy. In this 
expression, Ovib gives the per-particle free energy when 
the system is confined to an average basin of depth 0, 
and the basin enumeration function is defined such that 
C exp [Na{(j), p)] dip gives the number of inherent struc- 
tures with potential energy per particle <j) ± \d(f). For 
the supercooled liquid, both the vibrational free energy 
and basin enumeration function account only for those 
inherent structures and basins which are sufficiently de- 
void of crystallites. The partition function Z therefore 
corresponds to a restricted ensemble for the metastable 
liquid ll. 



The integral in Eq. |3.l| is performed over the per- 
particle potential energy of inherent structures and 
ranges from the lowest-lying amorphous configuration 
(with energy 0,nin) to the least stable arrangement (with 
energy 0max-) Presumably at these extreme energies, 
there exists a single, unique amorphous basin; it there- 
fore follows that cr ((/)min) — o'(</'max) = 0. The condition 
of zero configurational entropy does not exclude the exis- 
tence of inherent structures with sub-exponential system 
size dependence. Consequently, the notion of a unique 
amorphous inherent structure at cr = refers to the sub- 
set of basins that are thermodynamically relevant in the 
large-system limit. At this point, we should note that in 
the following discussion, the notation davihi(3, (j)* , p) /d4> 
is used to indicate the (/)-derivative of Cvib taken at con- 
stant values of its remaining natural variables, /3 and p, 
and evaluated at </>*. 

The transformation in Eq. offers enormous sim- 
plification in the analysis of liquids; given a and a^ih, 
for which good theoretical functionalities can be writ- 
ten, the usual SA'^-dimensional integral over positions 
in the configurational partition function becomes a one- 
dimensional integral over potential energy. Furthermore, 
Eq. 3.1 permits the separate determination of vibrational 
and basin contributions. Often the vibrational com- 
ponent can be successfully modeled using a harmonic 
approximation, while basin degeneracies can be under- 
stood through packing considerations jl^, ^ H^, H) . 
This casting of the partition function offers a particularly 
convenient framework for deeply supercooled liquids, in 
which intrabasin equilibration occurs on much shorter 
timescales than interbasin hopping . 

We focus our attention on the implications of the PES 
for supercooled liquids and ideal glasses. In the large 
system limit, the integral in will be dominated by 
a maximum value of the exponential term at cj) — cj)*. 
This means that, at the given temperature and den- 
sity, the system will sample configurations whose over- 
whelming majority have energy (/>*; thus one can con- 
sider only those contributions corresponding to a{(jf,p) 



and avib(/3, 

Z « ( 



,JV[ct(0* ,p)-/30* -/3avib(/3,0* ,p)] 



A/N « cp* -kTa{cp\p)+a,,^,{l3,<t>\p) (3.2) 

where there is a clear distinction between the basin 
((/)* — kTa) and vibrational (ovib) contributions to the 
free energy. In particular, Nka gives the configurational 
entropy, that is, the entropy due to the system's explo- 
ration of multiple basins of a given depth. This is con- 
sistent with conventional notions of entropy; when the 
system is confined to a single basin, as in a crystal, the 
configurational entropy vanishes ]40[| . 

The condition that determines the mean inherent- 
structure energy (f)* requires special attention. I n t he 
usual approach, the maximum of the integrand in 3T is 
located via the derivative of the exponential term 



da{4>\p) 



= /3 + /3 



9avib(/3, 0*,p) 



(3.3) 



i?(/3,0*,p)- 

This is an implicit expression for 0* as a function of (3 and 
p. Here we have introduced the function B for notational 
convenience; it appears frequently in this theory and is 
given by 



B[f},cP,p) = 



1 



9avib(/3, <^,p) 



(3.4) 



For simplicity in this analysis, one might assume the vi- 
brational free energy to be independent of inherent struc- 
ture energy, i.e., B((3, 4>, p) = 1. This is the rigorous limit 
of B at absolute zero, where the vibrational free energy 
and its (/)-derivative vanish. At finite temperatures, how- 
ever, this approximation gives meaningful results only 
when the curvature of basins is weakly dependent on their 
depth, a condition that might be expected for very deep 
basins (and low temperatures) where the majority of the 
particles are "well-packed." Still, we emphasize that this 
remains an approximation as a nonzero dependence of 
basin shape on depth has been found in a number of sys- 
tems [E6lEll0,0. 



The condition described by Eq. 3.3 holds for all tem- 
peratures if the slope of the basin enumeration func- 
tion tends toward infinity at ip-mm- With such behavior, 
the mean inherent structure energy decreases continu- 
ously to (pram as the temperature approaches zero (i.e., 
as (3 — > oo.) For basin enumeration functions with a 
finite slope at their minimum energy, the inherent struc- 
ture energy reaches ^min at a nonzero temperature. At 
lower temperatures, the value (jf which produces the in- 
tegrand maximum will always occur at the low-energy ex- 
tremum. This condition, in which the system is trapped 
in a unique amorphous minimum-energy configuration 
(or in one of a sub-extensive number thereof), is the ideal 
glass (IG) transition. The p* condition is then given by 



s.t. 



2£l _ P_ 

~ B{ 

0min(p) 



BW:kp) for/3<AG(p) (3 5) 



for/3>/3iG(p) 
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where I3ig{p) [= 1/^^ig(p)] is formally determined by 



da'{(t>min, p) 
dcj) 



Pig 



-B(/3lG,0min,P)' 



(3.6) 



In this case, a system in equilibrium at temperatures 
below Tig would be confined to the lowest-lying basin 
available at that density, with energy </)min- Accordingly, 
a = f or the entire range of temperatures T < Tiq. In 
Eq. 3.5, it is important to observe that with reasonable 
molecular interactions, the (j) derivative of the vibrational 
free energy should remain bounded as T ^ 0. Eq. 3.5 is 



the general case for basin enumeration functions. In the 
limit d<7((t>inin ) /d <P oo, we ha ve that Tig ^ and the 
expression in il reduces to |3.3| . 

The preceding equations provide a thermodynamic 
framework for an ideal glass {a = 0), motivated by the 
energy landscape formalism. The transition to this state 
at Tig is marked by the sudden confinement of the sys- 
tem to a unique, lowest-energy amorphous configuration, 
that is to say, the configurational entropy becomes zero. 
By our notation, we strictly identify Tig as an ideal glass 
transition, though previous studies have named the point 
at which the configurational entropy of the liquid be- 
comes equal to the configurational entropy of the crystal 
a Kauzmann point 0, ^ ^ Q, extending the 
original definition of "Kauzmann transition." The reader 
is reminded, however, that coincidence with a true Kauz- 
mann transition (i.e., as originally defined) occurs only 
when the vibrational entropies of the liquid and crystal 
are equal. To prevent any confusion, we avoid calling 
the ideal glass transition line a Kauzmann line. Further 
in this discussion we will explore the extent to which 
the observations regarding Kauzmann points presented 
in Sections || and || also hold for ideal glasses, which will 
effectively be a test of the crystal-liquid vibrational en- 
tropy equality. 



B. Equation of state 

It is apparent from the energy landscape formalism 
that there are two classes of basin enumeration functions 
p8[ . Those which have a finite slope at ipmin give rise to 
an ideal glass at finite temperature; those whose slope is 
infinite lack such a transition. To consider the implica- 
tions of these two possibilities on the equation of state, 
we examine the volume derivative of Eq. |3.2| . In the ideal 
glass-free case, this is 



(3 + f3 



i9avib (/?>*, p) da{(f)*,p) 



(3 dp 

2 ^Qvib(/3,0*,p) p^ da{(f)*,p) 
dp [3 dp 

PM0,p) + PisiP,p) 



(3.7) 



where the term in brackets vanishes according to the 
equilibrium condition in Eq. 3.3, and Pvib and Pis, 



the vibrational and inherent structure pressure contribu- 
tions, are conveniently defined by adding and subtracting 
{p^B/P)da/dp to the last two terms in this equation and 
rearranging to obtain 



PMP,p) 



Pis{P,p) 



, da^ib{f3,(j}*, p) 
dp 

[BiP,r,p)-l] 



da{(j)\p) 



--P(/3,<^ ,p)^^. 



dp 



(3.8) 
(3.9) 



This separation of the pressure components is con- 
structed so that the inherent structure pressure coincides 
with that measured in simulation studies of energy land- 
scapes. Perhaps more revealing, Pig is equivalently de- 
fined by the negative volume derivative of average inher- 
ent structure energy at constant configurational entropy. 
The equivalence with Eq. 3.£ is established by the math- 
ematical identity 



da_ 

dp 



da 



dcj)^ 
dp 



(3.10) 



when used in combination with the equilibrium condition 
in Eq. 3.3 and substituted in the pressure expressions in 



Eqs. 3.8 and 



Pvib(/3,p) = p' 



da^ihjP , cj)* , p) 
dp 

da,ibiP,r,p) (d_£_ 
dp 



PisiP.p) = p" 



dp 



(3.11) 
(3.12) 



where the 0*-derivative is taken along a path whose con- 
figurational entropy is consistent with the given density 
and temperature. For a more detailed discussion of the 
separation of vibrational and configurational contribu- 
tions to the pressure, the reader is referred to the ap- 
pendix. 

In the case of a finite-temperature ideal glass, Eq. 3/7 
still describes the pressure for T > Tig(p). Below the 
ideal glass transition, however, both (jf and a are con- 
stant at fixed density, constrained to their value at Tig. 
Using these constraints in the expression for the free en- 
ergy and subsequently taking its volume derivative then 
gives the pressure below the ideal glass transition: 



P 



Sflvib (/3, 0min , P) , Sflvib (/3, 0min , p) d(t>min (p) 



'P 



dp 

2rf0min(p) 



dp 



dp 

PU^{/3,p)+Pls{p) 



(3.13) 



where P^^^^ and P/g correspond to the first and last terms, 
respectively, and the prime symbols indicate the formal 
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diff erence in these definitions from those in Eqs. 3.8 and consider basin enumeration functions of the general form 



3.9 . Despite dissimilar expressions above and below Tiq, 
the individual pressure components remain continuous. 
Along the ideal glass transition locus, Tiq{p), the configu- 
rational entropy remains constant at zero and 0* — 0min- 
Considering this beha vior in t he alternate pressure ex- 
pressions in Eqs. 3.11 and 3.12, we have the equalities 



Pis(/3iG,p) = Pisip)- 



(3.14) 
(3.15) 



It is important to note here that the inherent structure 
pressure has no temperature dependence below the ideal 
glass transition, as is also the case with the inherent 
structure energy. Using the equivalence of pressure ex- 
pressions at Tig, one can write the inherent structure 
pressure for T < Tig in terms of its properties at the 
ideal glass transition: 



^l's(p) = -i^^(AG,0min,p) 
PIG 



d(7{(t)min,p) 

dp 



(3.16) 



With the equation of state in 3.13, we can immedi- 
ately evaluate the Sastry density, which we now define in 
a precise way as the zero-temperature minimum in the 
inherent-structure pressure as a function of density, 



Ps- 



(Ps) , r,d(pmin{Ps) 



dp 



(3.17) 



The implication of this expression is that Sastry den- 
sity is completely determined by the density dependence 
of the energy of the lowest-lying basin. This is not all 
that surprising, since one woul d ex pect A/N (f>min as 
T ^ and (j{(j)min, p) — 0. Eqs. 3.9 and 3.13 also provide 
the T > behavior of the pressure minimum and a cor- 
responding generalization to a finite-temperature Sastry 
density, which allows the correspondence with computer 
simulation studies. We note that the Sastry density de- 



fined by Eq. 3.17 also coincides with the density at which 
the inherent structure pressure along the ideal glass tran- 
sition locus has a minimum, and is consistent with the 
zero-temperature minimum of the total pressure, where 
any contribution from the vibrational free energy van- 
ishes. 



C. Basin enumeration functions 

To understand more specifically the implications of a 
liquid's PES, we now turn to a class of simple expressions 
for the basin enumeration function which have generated 
considerable interest [|4[ |2^, . Such generic func- 
tions must of course be taken as approximate if they are 
to be universally applicable. Nonetheless, these expres- 
sions provide a starting point for the analysis of super- 
cooled liquids and have even produced quantitatively ac- 
curate models in several simulation studies |l^, We 



cr(0,p) 



)(p)/N with u: 



a{p) 



(3.18) 



where f[u] is a dimensionless function varying between 
and 1, (Too is the maximum value of the basin enu- 
meration function, and u is the energy parameter made 
dimensionless by the inherent structure characteristic en- 
ergy range a and offset by an energy 0oo- We define the 
parameters 0oo and a such that u = — 1 for cj) = (pnim 
and (j)oc is the energy at the maximum value of the basin 



enumeration function (i.e., a 



J. For dis- 



tributions whose energy range is symmetric about (poo, 
we then have — 1 < u < 1 ((l>mhi < (Z* < '24>oo — (f'mhi)- 



Furthermore, the form in 3. IS gives the total number of 
basins from exp {Naoo) |19| . 

The above approach allows the "extraction" of generic 
properties of the basin enumeration function, such as 
height, width, and mean, from its specific functional 
form, given by f[u\. One of the simplest and most com- 
monly used models is the Gaussian landscape [p^ , given 
by 



/g 



aussian 



(3.19) 



As evident from this expression, the Gaussian form of 
the basin enumeration function has a finite slope at its 
energy minimum, and therefore, always leads to an ideal 
glass. 

For basin enumerations of the form in 3.1^ , equilibrium 
is determined by 



u*{j3, p) s.t. 



df[i, 



,u - B(/3.r,;)L(p) for/3</5iG(p) 
^u* = -l for (3 > Picip) 

(3.20) 

where u* is the equilibrium value of u and the glass tran- 
sition temperature is given by 



a{p)P] 



IG 



df[-l] _ 

-B(/3lG,?!>min,p)croo(p) du 



(3.21) 



When the vibrational free energy is independent of basin 
dept h {B = 1) , a scaling relationship evolves from Eqs. 
3.20 and [3.21 between the temperature, total number 
of basins, and basin distribution breadth. That is, the 
density dependence of the mean inherent structure en- 
ergy is contained in a dimensionless temperature vari- 
able a(3/ (Too- Galculati ng t he in heren t structure pressure 
contribution from Eqs. 1^ and |^ 



Pis(/3,p) - p' 



-Xp) I ^* da{p) 
dp 



dp 

p^BiP,r,p) 



, dcFoo [P] 



dp 



(3.22) 



which holds above the glass transition. As T ^ Tig, 
the last group of terms in this equation vanishes since 
f[u*] — » 0, and the first group on the right hand side 
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becomes the same as that in Eq. 3.13 by the definition of 
6r^ and a. 



The form in Eq. [3.1^ has interesting imphcations for 
the apparent coincidence of the Sastry density and the 
zero-temperature Kauzmann point when the crystal and 
liquid vibrational entropies are equal. We assume first 
that we have a functional form for the basin enumera- 
tion function which gives rise to an ideal glass transition, 
i.e., df[—l\/du is finite. Then using the fact that in the 
T = limit the vibrational free energy is independent of 
basin depth and following Eq. 3.21, the Tic-ps conver- 
gence implies 



a(ps) (df[-l] 



du 







(3.23) 



where ps is found from Eq. 3.17. This result implies 
that either cr oo diverges or a vanishes at the Sastry den- 
sity. The first case appears unlikely; elementary argu- 
ments about the shredding behavior of low-density inher- 
ent structures imply that CTqo should vary continuously 
across ps . It thus appears that the range of inher- 
ent structure energies, given by a, must shrink to zero 
at the Sastry density. Qualitatively, this seems plausible. 
At smaller volumes, the system boundaries can serve to 
stabilize poorly packed, high-energy inherent structures, 
but as the density is lowered into a state of tension, this 
is no longer the case and the energy of the least stable 
structure necessarily decreases. With the increasingly 
constraining requirement of homogeneity, the maximum 
attainable inherent structure energy may be forced to 
converge on the minimum as the Sastry density is ap- 
proached, that is, 0,„ax 4>min at PS. 



IV. ELEMENTARY LANDSCAPE MODEL 

We now derive properties of a basic energy landscape 
for simple liquids which serves as a minimum descrip- 
tion for their analysis in the landscape formalism. We 
consider a simple model which is particulary convenient 
to analyze in the energy landscape paradigm. This ide- 
alized system consists of structureless, spherically sym- 
metric particles interacting in pairs through soft-sphere 
repulsive forces of the type r~" and each experiencing 

2^, 0. The 
es at a given 



a density-dependent mean-field attraction [ 
potential energy of a configuration of partic 
density is therefore given by 



U{r^^,p)^\e{a/n,r-Nap 



i<j 



(4.1) 



where r"^ represents the positions of the particles, e and 
a are the characteristic energy and length scales of the 
soft-sphere part, n is the soft-sphere exponent which con- 
trols the degree of the repulsion, rij is the distance be- 
tween particles i and j, and a, a positive number, is the 
mean- field parameter with units of [energy] x [volume]. 
There are in fact only two independent parameters in this 



expression, the soft-sphere and mean-field coefficients, 
ecr" and a. This particular choice of the potential en- 
ergy function is motivated by the observation that liq- 
uid structure is primarily determined by repulsive forces 
while attractive interactions can be successfully incorpo- 
rated by a background potential field serving to hold the 
particles together Q . The advantage of this system 
is that the repulsive part of the potential energy in any 
configuration is determined by the dimensionless particle 
positions, s — V~^^^r, in such a way that the volume 
dependence can be readily extracted: 



U/N = e(7"y-"/3iV-i^s-" 



ap 



7(s3^)p"/3 _ ap 



(4.2) 



where 7(5'^^), with units of [energy] x [volume]"/^, has 
been defined to give the portion of the soft-sphere term 
independent of volume, and is given by 



l<j 



(4.3) 



Due to this kind of volume dependence, the nonideal soft- 
sphere part of the configurational partition function de- 
pends on a single scaling variable z = {(3e)^^^pa^ [^ . 
This reduction in parameter space has made the soft- 
sphere/mean- field (SSMF) model attractive for studies 
of the liquid state [|2^, |3^ . Previous investigations using 
an equation of state based on reasonably accurate sim- 
ulation data have shown in the model the presence of a 
Kauzmann locus and the intersection of this curve with 
the T = liquid spinodal Q. We aim to extend the 
work in these studies by using the scaling embodied in 
Eq. 4.2 to develop a model energy landscape. 



It is well known that the total number of inherent 
structures in the family of soft-sphere models is indepen- 
dent of density |3^, This follows directly from the 
scaling effect of the volume on the potential energy and 
forces, which means that inherent structures are defined 
by steepest descent trajectories in the dimensionless par- 
ticle positions, s^^ . Extending this observation, the ener- 
gies of the most and least stable inherent structure config- 
urations must scale with density as p"/"^. Therefore, the 
number of inherent structures within any fractional por- 
tion of this energy range must be density- independent. 
The result is that the functional form of the basin enu- 
meration function for soft-sphe res is rigorously density- 
independent, that is to say, Eq. 3.18 is exact in this case. 
(We should note here that these observations are also the 
limiting high-density case for more realistic pair potential 
functions which incorporate an inverse-power repulsion, 
such as the Lennard- Jones potential. In such systems 
at high densities, the relevant configurations of the par- 
tition function sample only the repulsive part of these 
potentials.) Finally, the appended mean-field term does 
not change the inherent structures or their number, but 
simply serves to shift their energy. With these consider- 
ations, we can elaborate on some of the properties of the 
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FIG. 2: Sastry density for the soft-sphere plus mean-field 
model as a function of repulsive exponent, n. The limiting 
case for n — + <x, shown by the dashed line, is ps/po ^ 1- 



basin enumeration function: 



'min(p) 




(4.4) 






(4.5) 


a{p) 


= (Too - 7min)p"/^ 


(4.6) 


O-oo(p) 


= Coo 


(4.7) 



where 7min and 700 correspond to the values of 7 for con- 
figurations at the low energy extreme and maximum of 
the soft sphere basin enume ration function, respectively. 
Using expressions 4.4-4.7, the Sastr y den sity for the 



SSMF system follows directly from Eq. 3.17 



PS 



1 /n\ /n 
2 \i) V3 



(4.8) 



where we have introduced a characteristic density poi de- 
pendent on the soft-sphere exponent n and defined by 



Po 



Tniir 



(4.9) 



The significance of po stems from Eq. 4.8; in the limit of 
hard sphere repulsive interactions for which n — > 00, the 
Sastry density tends towards this reference density. This 
trend is depicted in Fig. which shows the evolution of 
PS over several orders of magnitude in n, starting with 
n = 4. Despite the existence of a zero-temperature pres- 
sure minimum, one important distinction between this 
model and ones incorporating more realistic attractive 
interactions is the absence of a distinguishing physical 
change at the Sastry density. The SSMF system remains 
completely homogeneous below ps a-nd does not exhibit 



the "shredding" behavior normally associated with this 
regime. This makes it difficult in this case to justify, 
at least in a qualitative sense, the notion that the Sastry 
density represents a point of zero configurational entropy 
due to the constraint of homogeneity ||2^ . Unfortunately, 
this type of justification requires a model of greater com- 
plexity, for which there is currently no obvious candidate. 

In order to evaluate the equilibrium pressure in the 
SSMF fluid, we model the vibrational free energy in the 
harmonic approximation. For consistency, we use the 
classical form @ |4|, ^ |§: 



avib(/3,'^,p) 




(4.10) 



(4.11) 



where is taken as an Einstein temperature, h is 
Planck's constant, m is the particle mass, and the term in 
brackets is the geometric average of the eigenvalues of the 
Hessian matrix (i.e., the matrix of second derivatives of 
the potential energy with respect to particle coordinates.) 
In this expression, there is an implicit dependence of 
on basin depth and density, both whose effects manifest 
in the second derivative of the potential. This means 
that for given values of (/) and p, the Einstein tempera- 
ture is essentially derived from a "representative" inher- 
ent structure of that basin depth at the specified density. 
Making a transformation similar to that in Eq. 4.2, we 
have 



Qr = 




p("+2)/3 



7Vl/3^\ 
/ 



p(»+2)/3 



(4.12) 



where a function of the scaled basin depth u, gives the 
density-independent part of the second derivative term. 
Essentially g contains any dependence of basin curvature, 
and hence vibrational free energy, on basin depth; a (fy 
indcpendent Ovib simply implies that g is constant. Using 
the harmonic expression for the vibrational free energy 
in Eq. 3.20, the ideal glass locus is given by 



Ag(p) 



1 



'~ioo Tmin 



df[-i] 

du 



3dln5[-l] 
2 du, 



Cp- 



-n/3 



(4.13) 



where we have condensed the constants in this expres- 
sion into a single constant C for clarity. The equivalent 
expression for temperature is 



Tig(p) 



1 



n/3 



(4.14) 



The surprising consequence of this result is that an ab- 
solute zero ideal glass transition is reached only in the 
limit of infinitesimal density. This means that the Sas- 



try density, which is finite via Eq. 4.8, does not intersect 
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the ideal glass curve at T = 0. This result may seem at 
odds with the results of Ref. in which it was found 
from numerical equations of state for the n — 9 SSMF 
model that the zero-temperature limit of the Kauzmann 
curve intersected the terminus of the liquid spinodal. One 
must remember, however, that the focus of Ref. ||2^ 
was the Kauzmann locus whereas we strictly consider 
ideal glasses; the difference in the approaches is of course 
whether the total or configurational entropies are of con- 
cern. Based on the discrepancy between the two studies, 
apparently there are subtle but fundamental differences 
in the implications the definitions of these two phenom- 
ena have for the analysis of liquids. The present work 
therefore emphasizes that an ideal glass transition does 
not necessarily correspond to a Kauzmann point; that is, 
the sudden confinement of a liquid to a unique, lowest- 
lying amorphous basin in its potential energy landscape 
can be a distinct occurrence from S'uquid = '5'ciystai- 

Turning now to the pressure along t he id eal g lass tran- 
sition locus, insertion of Eqs. 4.4-4.7 and 4.1C into Eq. 
3.13 yields 




1.00 



^is(p) 



(4.15) 



for the inherent structure contribution, and 



■Pvib(/3lG,P) = 



n + 2 



2 /3ig(p) 
n + 2 



2C 



-P 



i/3+l 



(4.16) 



for the vibrational part. Without the functional form for 
/[it], the expression for the pressure above Tig is unspec- 
ified in the model. Though thorough discussion of such 
specific expressions is beyond the scope of this paper, we 
present for the purposes of demonstration the equation 
of state formulated according to a Gaussian landscape 
(see Eq. B.19). To give the most basic expression, we 
make the assumption of constant g, i.e., that the vibra- 
tional free energ y is independent of basin depth. It then 
follows from Eq. 3.7 that 



-ap 



loo - l22fL_ZlIli!l)!/3p"/3 

2.(Toa 

n + 2 p 



2 (3- 



(4.17) 



From this equation, one can see more clearly the possi- 
bility of a T-independent inherent structure equation of 
state. The inherent structure pressure, given by all but 
the last term in 4.17, has a weak temperature dependence 
if the quantity cTooT is reasonably large. In this case, the 
term in brackets is essentially given by the constant value 
7oo and the inherent structure pressure approaches 



IS 



370 



2/3+1 



ap 



(4.18) 



Though the SSMF system is missing any dependence 
of the total number of basins on density, which would 



FIG. 3: Pressure isotherms for the landscape-derived equa- 
tion of state of the soft-sphere/mean-field system, fitted to 
Lennard- Jones state points. The predictions of the theory 
(lines) correspond to T = 0.7 — 1.3 in increments of 0.1, from 
bottom to top. The filled diamonds are the simulation results 
from and the squares are the corresponding predictions 
of the theory. All variables are expressed in reduced Lennard- 
Jones units. 



almost certainly accompany explicit attractive interac- 
tions, one could imagine this to have a relatively small 
effect at commonly-studied temperatures. Qualitatively, 
then, the inherent structure pressure would continue to 
appear temperature-independent at high T . 



We have fitted the equation of state in 4.17 to the sim- 
ulation results for the Lennard- Jones system in Ref. [^ . 
In Eq. 4.17, the parameters 7min and Cqo are not inde- 
pendent because they appear in the same coefficient. In 
fitting the model, we set n = 12 and minimize the error in 
predictions for the pressure at several of the high-density, 
low-temperature state points in |Q, obtaining values of 
4.64, 16.5, and 0.763 for 700, a, and (700 - 7min)^/o-oo, 
respectively. Ou r results are displayed in Fig. |^. In 
general, Eq. 4.17 predicts the pressures well at increased 
densities, which is the expected result since at these state 
points, the Lennard- Jones system samples primarily the 
repulsive part of its potential. Interestingly, though, the 
predictions fare relatively well for the two highest densi- 
ties even at temperatures near the critical temperature. 
This is somewhat surprising considering the approxima- 
tions in the theory, including the harmonic approxima- 
tion and the assumption of a b asin-independent vibra- 
tional free energy. Using Eq. 4.8, we also obtain a Sastry 



density of 0.76 for this system; though this value differs 
from the simulation result of 0.89 in Ref. ||l^, one would 
expect a deviation since the SSMF model neglects any 
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explicit attractive interactions which are sure to be im- 
portant at pg. 



V. CONCLUSIONS 

We have presented a theoretical framework for liquids 
which incorporates the statistical properties of an energy 
landscape, is capable of describing an ideal glass transi- 
tion, and corresponds with simulation studies performed 
in the inherent structure formalism. In particular, the 
theory provides an explicit expression for the equation 
of state of a liquid and formally separates the pressure 
into vibrational and inherent structure components, both 
above and below the ideal glass transition. Using the 
pressure separation, we have shown the presence of the 
Sastry density and its connection to the liquid spinodal. 
Finally, we have used the theory to develop an elemen- 
tary model of an energy landscape based on soft-sphere 
particles interacting with an additional mean-field attrac- 
tion. The model appears reasonably descriptive of simple 
liquids, suggesting the usefulness of a landscape-based re- 
formulation of supercooled liquid thermodynamics. An 
important conclusion from this work is that the ideal 
glass and Kauzmann loci are quite distinct for the soft- 
sphere model with mean-field attraction. 

Though the current work provides a detailed theoreti- 
cal picture for energy landscapes, a number of questions 
remain. Improved knowledge of the vibrational free en- 
ergy at low temperatures is needed in order to provide 
a microscopic understanding of the distinction between 
an ideal glass transition and a Kauzmann equal-entropy 
point. The apparent coincidence of the liquid spinodal 
with the Kauzmann locus remains unresolved; the re- 
sults of our theory suggest that notions of an ideal glass 
at zero configurational entropy fail to capture what is oc- 
curring at this presumed convergence. Finally, because 
our theory is based on the energy landscape, it contains 
information useful for both thermodynamics and kinetic 
phenomena (e.g., diffusion rates) j^, |2^. A useful fu- 
ture investigation might examine its ability to capture 
the slowdown of kinetic processes near the glass transi- 
tion. 
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APPENDIX: IDENTIFICATION OF THE 
INHERENT STRUCTURE PRESSURE 

The par ticul ar s eparation of the pressure components 
in Eqs. 3.8 and 3.9 is constructed to coincide with simula- 



tion studies of energy landscapes, in which inherent struc- 
tures are found by energy minimization to local poten- 
tial energy minima and their pressure determined from 
the usual virial expression. The algorithmic procedure in 
these simulations may be thought of as an instantaneous 
quench to a non-equilibrium state at zero temperature. 
Using a recently proposed thermodynamic formalism for 
such systems |39|, we define the inherent structure pres- 
sure as the negative volume derivative of the free energy 
of the quenched state. The non-equilibrium, per-particle 
free energy is given by 



aNE(/3int,/^q,p) = - /3i„tcr(0*,p) -I- avib(/3q,0*,p) 

(A.l) 

where /3int refers to an internal temperature character- 
izing the out-of-equilibrium dynamics between inherent 
structures, /3q refers to the temperature to which the sys- 
tem is quenched, and (jj* is the equilibrium value of (f) 
prior to the quench, a function of the original temper- 
ature of the system p. In the case where the quench 
temperature is absolute zero, the expression for Pi^t is 



Pint — 



B{p,4>\py 



(A.2) 



Using this expression in the volume derivative of Eq. A.l 
and taking the limit T — > leads directly to the ex- 
pression for inherent structure pressure in Eq. 3.9. The 



vibrational component of the pressure then follows from 
the remainder, P — Pig. 

An alternative but equivalent definition of the inherent 
structure pressure is the negative volume derivative of in- 
herent struct ure e nergy at fixed configurational entropy, 
given by Eq. ^.12 , In the case of soft-spheres interacting 
through an inv erse-power potential, one can demonstrate 
that Eq. [3.12 is the correct Pis, that is, the one actually 
measured in simulations. The unique property of soft- 
spheres is that upon isotropic compression in which all 
particle coordinates are simply rescaled, a configuration 
that is a potential energy minimum remains a minimum. 
As a result, the number of members in a given ensem- 
ble of such configurations does not change with this type 
of volume deformation, and hence the same is true for 
the configurational entropy corresponding to any inher- 
ent structure. Since the virial measures the response of 
the potential energy to an isotropic volume change of 
this sort and because of the constancy of configurational 
entropy during the compression, simulation pressure cal- 
culations for soft-sphere inherent structures are indeed 



given by Eq. 3.12 and hence also 3.9 
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